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ABSTRACT 

Encouraged by the recent construction of fuzzy sphere solutions in the ABJM theory, we re-analyze 
the latter from the perspective of a Matrix-like model. In particular, we argue that a vortex solution 
exhibits properties of a supergraviton, while a kink represents a 2-brane. Other solutions are also 
consistent with the Matrix-type interpretation. We study vortex scattering and compare with graviton 
scattering in the massive ABJM background, however our results are inconclusive. We speculate on 
how to extend our results to construct a Matrix theory of ABJM. 
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1 Introduction 

Based on the earlier Bagger, Lambert & Gustavsson (BLG) construction [1, 2, 3, 4] of an M2-brane 
action in terms of 3-algebras, Aharony, Bergman, Jafferis k, Maldacena (ABJM) [5] constructed an 

M = 6 supersymmetric action for the IR of N M2-brancs probing a C^/Zk singularity, a U{N) x U{N), 
level-A; Chern-Simons-matter gauge theory whose fields transform in the bifundamentals. Not only 
was this the first time that an action for an arbitrary number of M2-branes was written down but it 
also allowed further insight into the structure of M-theory, as well as a new and exciting example of 
the AdS/CFT with its gravity dual being the large k Umit of AdS4^ x S'^/Zk, i.e. AdS4 x CP^. 

But while most of the work on the ABJM model concentrated either on the field theory side, or on 
the AdS / CFT duality, it is worth bearing in mind that the original interest in the multiple M2-brane 
system was, of course, the potential for a better understanding of M-theory. Indeed, if we are to use 
the ABJM model toward this end, the most natural possibility that comes to mind - since we are, 
afterall, dealing with a gauge theory oi N x N matrices - is a Matrix theory-type construction. In 
fact, very much in the spirit of M-theory, it was found in [6, 7] that the BPS fuzzy funnel solution 
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of pure ABJM, or the vacuum fuzzy sphere solution of the massive deformation of ABJM gives rise 
to a D4— brane on 5^, together with the correct small fluctuations action, in the classical limit. This 
D4— brane appears as an M5-brane on /Z^ at large /c, in a similar way to how D2— and -D4— branes 
appear in the Matrix theory of Banks et.al. (BFSS) [8, 9, 10, 11]. 

It seems reasonable therefore to expect that a Matrix theory-type model can be constructed 
out of the ABJM model. Following the logic of Matrix theory, we need to find a classical solution 
of ABJM corresponding to a spacetime supergraviton. Since such a solution must be localized on 
the worldvolume, as well as in transverse space, it must be a particular type of vortex solution. 
However, since the ABJM action, unlike its BFSS counterpart, is conformal, we will see that a better 
definition of its Matrix model is given by a maximally supersymmetric deformation of ABJM [12]. 
Consequently, after an initial analysis of pure ABJM model, for most of the rest of the paper we 
focus on its maximally supersymmetric deformation. After identifying various spacetime branes as 
classical solutions of the model, we compute supergraviton scattering and compare this with vortex 
scattering in ABJM. Unfortunately, we find a mismatch between the simplest calculations on both 
sides that manifests in the associated interaction potentials. We will argue that this mismatch arises 
because a naive application of the BFSS model does not identify the correct calculations on both 
sides that are supposed to match. We will then speculate on the Matrix theory rules for the correct 
identification of the two sides. 

The paper is organized as follows. In section 2, in the interests of being self-contained, we provide 
a lightning review of BFSS theory, focusing on those elements that carry over to our case, and then in 
section 3 we will describe in detail our set-up and what we expect to find. Section 4 is devoted to an 
analysis of solutions of pure ABJM theory, and we argue that, while it is possible to identify a vortex 
solution with a supergraviton, such a solution is practically unfeasible due in no small part to the 
infinite energy of the corresponding background D2-brane. Instead in section 5 we will focus on the 
massive deformation of ABJM and identify its known solitonic spectrum (consisting of brane-filling, 
kink and vortex types) with branes in spacetime. In the process we identify the background in which 
the M2-branes corresponding to the massive ABJM move. In section 6 we calculate the scattering of 
supergravitons in this massive background, and in section 7 the corresponding scattering of vortices 
in massive ABJM. Finally, after noting a mismatch of the calculations of sections 6 and 7, in section 
8 we speculate on the possible definition of the sought-for ABJM Matrix theory. 



2 A short review of the BFSS matrix theory 

The BFSS Matrix theory of [8] is based on a discrete light cone quantization (DLCQ) of M-theory 
with a compact circle. In the heuristic derivation given by Sen and Seiberg [13, 14], this light- 
like compactification of the given M-theory is related to a space-like compactification of a different 
M-theory in a decoupling limit, in which the only thing that remains is a decoupled theory of N DO- 
branes. Each DO-brane in this description corresponds to a single unit of momentum in the compact 
(11th) M-theory direction. Its action is given by 
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where the = X^/gl'^ are nine scalars corresponding to the nine transverse directions, = dt+iA, 
and ^ is a 0+1 dimensional U{N) gauge field. In addition R is the compactification radius, Mp, 
the 11-dimensional Planck mass and the 9 are fermionic superpartners of the that transform as 
spinors under the 5*0(9) group of transverse rotations. 

The simplest classical solution of the BFSS Lagrangian 

X\t) = {xi + vH)lMxM, (2.2) 

can be understood as M < DO-branes located at xq and moving with velocity v. In spacetime, 
this corresponds to a pointlike object with M units of momentum along the 11th direction^ and is 
interpreted as a supergraviton. Other classical solutions correspond to other D— branes in spacetime 
- D2— branes and D4— branes - with different possible geometries. 

Arguably, the calculation that received most attention in Matrix theory, ostensibly providing the 
first real test of the BFSS conjecture, is the matching of the interaction potential of two spacetime 
supergravitons with the corresponding interaction potential for two corresponding objects of the form 
(2.2) in Matrix theory. The calculation of the supergraviton interaction potential, as described in 
[15], is based on the observation of 't Hooft [16] that the tree level Rutherford scattering interaction 
potential in gravity (mediated by single graviton exchange) can be calculated by scattering two 
gravitational Shockwaves. One of these is described by an Aichelburg-Sexl shockwave - which serves 
as a heavy source graviton, while the second is a plane wave probing it. Note that the 't Hooft 
calculation was in a flat space background, but this procedure was also applied successfully to the 
curved space case [17, 18]. 

The first step in determining the interaction potential in the more general case, then, is to calculate 
the Shockwave in the given spacetime background. This plays the role of the graviton wavefunction 

and corresponds to adding an h (dx~)^ term to the background metric. These 'pp' Shockwaves 

have the remarkable property (not shared by many solutions in the highly nonlinear field equations 

of general relativity) that the linearized solution is exact, i.e. h is the solution of the Poisson 

equation 

Ai-g^{x)h..{x) = Qd^{x), (2.3) 

where the source Q depends on the momentum of the wave in the 11th dimension, pu = Ng/R and 
Mp. In the BFSS case, there are nine transverse directions (and two parallel directions - time and 
the 11th dimension, in which the wave propagates), so, by dimensional analysis, the solution must 

scale like 1/r^. More precisely h = ^^ 5{x~), and the source graviton is in a state of definite 

P- so that averaging over x~ G (0, 2'kR) gives 

h -J^l (24) 

To compute the interaction potential we then scatter a probe scalar (plane wave) off the shockwave 
background. Practically speaking, this requires considering the h-- part of the metric as an inter- 
action and using standard quantum mechanical perturbation theory. Consequently, the free wave 



^Or, equivalently, carrying DO-brane charge M. 
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equation = becomes 



h. 



-{do - m 



= V(f) (2.5) 



2 

where (p = 0e*^^" and p = im. At this stage, a few points deserve some elaboration. 

• Momentum in the 11*'* dimension acts hke a mass from the 10-dimensional point of view. 

• Even though gravitons move at the speed of light, we call a graviton moving entirely in the 11*'* 
dimension a f = solution. This will be an h gravitational perturbation. 

• The probe graviton has nonvanishing (i.e., a relative velocity between the probe and the 
source), and so will also propagate ever so slightly in the other ten dimensions. 

• For both gravitons, however, x" plays the role of time which is as it should be seeing as how 
we are in a DLCQ description. 

With this prescription (and the usual relativistic normalization), the probe graviton wavefunction 

where E"^ = fp' + m? , E = ^^^^ and p = J^^ - The 1-loop interaction potential in momentum 
space is computed through the S-matrix 

2T,5{Ein-Eout)V^^} = = j d^^'y 4>out{y)*V {y)4>in{y) ■ 
On substituting for the graviton wavefunction and the expression for V{y) above, 

^(1) ^ ^■(^"' ~"')" ^ f dxo e*(^-*-^-)^o / (fx e^(p---P-*)*/i__(f) , (2.7) 
SO that after integrating, the interaction potential can be read off as 



y«(x) = /.__i^^:.rn^/.__ + (9(.«). 
Finally, using m = Np/R gives the full one-loop result 



(2.8) 



Vintix) ~ NsNp-^ . (2.9) 

It is this interaction potential - in particular the coefficient of the v'^/r'^ term - that is then matched to 
the one-loop Matrix-theory of BFSS [15]. To compute the interaction potential of two supergravitons 

X\t) = {x{ + v\t)lN,xm + (4 + vit)lN,xN2, (2.10) 

separated by r = j.-ri — X2\ and with relative velocity v = \vi—V2\ in Matrix theory, we note that since 
the two diagonal blocks are non- overlapping these can be interpreted as two distinct supergravitons 
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with correspondingly different extent in tfie transverse directions. If we choose v\ = —v^ and = +v'^, 
we get that 

H = ^T.m-=^^^^±^'^ = ^ + ^ = E,+E,, (2.11) 

2i?'' R 2 2 2 1^^' \ J 

from the Matrix action. This is the free nonrelativistic energy of the supergravitons. To find the inter- 
action potential, we must go to one-loop, in which case one takes the one-loop fluctuation determinant 
around 

X\t) = {b' + vH)lNxN, (2.12) 

where now ¥ = x\ — and v'^ = v\ — ■ul are the relative positions and velocities respectively. The 
calculation of one-loop determinants amounts to just the zero point fluctuation X]„<jJn/2- Actually, 
this is the same calculation that one docs to compute the quantum mass of solitons like, for example, 
a kink. There one calculates this sum in the background of the kink solution (see, e.g., [19]). The 
result is in perfect agreement with (2.9). 

To summarize then, the leading order calculation of the interaction potential on the Matrix theory 
side matches the leading order result on the gravity side, even though the former is at one-loop, while 
the latter is classical. At this point, it is only natural to ask: 

Is there an analogous computation of interaction potentials that can he performed in ABJM? 



3 Set-up and expectations 

We saw that the BFSS model - a (O-l-l)-dimensional U{N) Matrix model on the worldvolume of N 
DO-branes - describes M-theory in discrete light cone quantization. However, since the Z?0— branes are 
momentum modes on the compact 11th dimension, this description of M-theory is not a fundamental 
one. Instead, as shown by Sen and Seiberg in [13, 14], it appears because of the equivalence of the 
original M-theory with a decoupled theory of DO-branes living in another M-theory. Any fundamental 
description of M-theory must involve M2-branes instead, but we don't know how to formulate it. In 
the large N limit classical I?— brancs appear as solutions of the classical BFSS theory. For example, 
D2-brane solutions in [8] were found, wrapping a fuzzy torus defined through, 

= RgP', = RgQ 
P = ^/Np; Q = ^fNq 

[Q,P] = 2iTi (3.1) 

where U = e'*' and V = e*^ are the "clock" and "shift" operators of the fuzzy torus satisfying 

UV = e^VU 

= V^ = 1 (3.2) 

Another instructive example is Matrix theory in a pp-wave background {i.e. the BMN Matrix model 
[20]). Here, 2— branes wrapping a fuzzy S"^ are also a solution. There are two ways to think about 
this: either as another example of the same BFSS construction as BFSS, only in a different spacetime 
background, or as a massive deformation of the BFSS model since the pp-wave corresponds to a 
mass deformation on the brane worldvolume. The presence of the mass deformation also serves to 
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better define the Matrix theory as it makes states discrete instead of continous. If Matrix theory is to 
correctly describe M-theory (and its dimensional reduction to type IIA string theory) then it should 
be able to describe all D— branes in the theory and not just D2— branes. For example, a Z)4— brane 
wrapping an S'^ was found in [9], following the earlier works of [10, 11], but the solution is not without 
several unresolved subtleties. In general, finding the complete spectrum of D— branes from Matrix 
theory remains a very difficult problem. The D2— and D4— branes already found are reductions to 
ten dimensions of M2— and M5— branes, and while they are a minimum necessary for the spectrum 
of M-theory, they are by no means sufficient. Indeed, we would also need to find a -D6— brane, coming 
from an eleven dimensional KK monopole, and a Z)8— brane^. 

Fortunately, the recent construction by Aharony et.al. [5] of an IR action for N M2— branes at a 
C^/Zfc singularity offers some much needed hope. Since this ABJM model is also a theory N x N 
matrices^ , it is natural to ask whether one can construct a new kind of Matrix theory from the ABJM 
model, perhaps one whose formulation is more fundamental instead of the somewhat derived one of 
BFSS. In support of this idea, it was shown in [6, 7] that D4— branes wrapping a fuzzy appear 
as solutions of (a massive deformation of) the ABJM model, in much the same way as in the BFSS 
Matrix theory. In as much as the BFSS Matrix theory is related to the maximally super symmetric 
BMN Matrix model by a massive deformation of the former, the pure ABJM model can be mass- 
deformed to yeild a maximally supersymmetric massive ABJM model [22]. In this case the fuzzy 
funnel solution of the pure theory stabilizes to a fuzzy sphere in the massive one. This hints then, 
that the massive deformation gives a better definition of a proposed Matrix model. Nevertheless we 
will first begin with an analysis of the pure ABJM model and then show how we are driven to its 
massive deformation. Our strategy will be as follows: 

• The presence of the D4 wrapping a fuzzy solution suggests that we start with a similar 
set-up to BFSS. 

• The first issue to be checked is whether there is a classical solution corresponding to a super- 
graviton, for which we can then compute scattering. 

• Thereafter, we will test whether other D-branes appear as classical solutions^. 

We immediately notice a difference from BFSS; namely that the N M2— branes of ABJM have two 
spatial worldvolume directions while a supergraviton needs to be pointlike in both the transverse and 
parallel space directions. It also carries DO— brane charge Af < so we must look for a vortcx-typc 
solution of ABJM, pointlike in the transverse space, and with some DO— brane charge AI. Once we 
find such objects, we can then try to scatter them, and match against scattering of supergravitons. 

Even having identified which objects to scatter, the issue of matching superpotentials is still a 
little murky since there is no a priori guide which term should correspond to which. This is not 
unique to our ABJM computations either. Indeed, in the BFSS model the leading supergraviton 

scattering is a tree-level (classical) interaction, while on the Matrix side the leading term is a one- 
loop interaction with the classical interaction vanishing. Moreover, on the gravity side the tree-level 

^The latter would appear in the massive type IIA string theory, for which the Matrix theory was constructed in [21]. 
Its properties remain largely unexplored. 

^Albeit one with a U{N) x U{N) gauge group and bifundamental degrees of freedom. 

''And reserve the right to not be surprised if, as for BFSS, it turns out to be difficult to find all of them. 
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interactions yield an infinite series that should match various higher loop correction terms in the 
Matrix side but the precise matching cannot be guessed before doing at least a computation of the 
Ns,Np and Mp dependences. 

In our case the increased computational complexity means that even the Ns,Np dependence is 
hard to obtain in the ABJM side while the Mp dependence cannot be guessed before computing the r 
dependence, as we discuss below. Our goal therefore will be to first compute and match the simplest, 
leading terms on both sides. On the gravity side, this will be a similar tree-level term as in the BFSS 
case, but now in the Matrix model, the leading classical interaction (for the vortex scattering) is no 
longer vanishing. 

One other difference between the ABJM and BFSS models is that pure ABJM is a conformal 

theory while the BFSS Lagrangian, having Mp as coupling, is not. In the BFSS interaction potential, 
the Mp dependence appears from two sources, one being as a coupling dependence, and the other 
as giving the unit of length when translating from BFSS to gravity. The latter is of the same kind 
as in AdS/CFT where the Ig dependence of gravity calculations also appears by introducing a unit 
of length in the conformal calculations of the field theory. Since pure ABJM is also conformal, only 
the latter type of Mp dependence would be available. So the matching is a priori more constrained 
than in the BFSS case. Of course, as we said above, the massive ABJM will be found to be more 
useful for defining a Matrix model. There we will see that we have one more parameter - the mass 
deformation // - potentially improving the situation. We will return to these issues in some detail in 
section 8 where we describe the systematics of matching in BFSS and how they might apply to our 
case. 

While on the topic of expectations, it is worth asking at this point what other D— branes we expect 
to find as classical solutions of the ABJM model? Certainly, as in BFSS, we would want at least an 
example of a Z)2— brane and a Z)4— brane, which would mean that M2— and M5— branes appear in 
the theory. These are already present however; for the D2— branes there is at least the configuration 
of the ABJM 2— branes themselves^ while for Z)4— branes there is (at least) the solution corresponding 
to a £)4— brane on a fuzzy S*^. Of course, this is by no means satisfactory, in part because we now 
have the peculiar fact that directions parallel and transverse to the ABJM worldvolume are different 
so we would like to see at least one example of a "shifted" brane. For example, we would expect 
to find kink solutions having a fuzzy one- dimensional transverse space and a fuzzy three-dimensional 
transverse space, corresponding to D2— and 154— branes respectively, each with only one direction 
parallel while transverse D2— and Z)4— branes with no parallel directions should be described as 
vortex solutions with fuzzy two-dimensional and a fuzzy four-dimensional transverse space. We will 
seek to find at least one example among these solutions. 

4 Pure ABJM solutions 

Solutions of the pure ABJM theory and its massive deformation appear to be related, usually in a fairly 
nontrivial way. For example in [6, 7], it was shown that not only does the maximally supersymmetric 
fuzzy sphere ground state of massive ABJM becomes the ^-BPS fuzzy funnel solution of pure ABJM, 

^See, for instance [23] for an example of how the worldvolume theory of N Z)2— branes arises. 
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but they also enjoy many shared properties, like the same unrescaled bosonic action for fluctuations. 
Certainly then, we would not be surprised if the vortex and kink solitons found in both pure and 
massive ABJM [24, 25, 26], turn out to be similarly related. In this section we will discuss the solitons 
of pure ABJM, though for technical reasons we will be forced to switch to massive ABJM in the next 
section. 

To this end, we begin by asking if we can find vortex solutions of pure ABJM theory that can be 
interpreted as supergravitons in spacetime? A good place to start answering this question is with the 
ansatz presented in [25] for a general class of vortex solutions, 

= vH, for A = 2, 3, 4, 

N-l 



a=l 



VM 



G{z) 



ddhi\ya\ 



'12 



F]o = — 



nN-l ' 
a=l Va 

)\y\yI 



(4.1) 
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where Y^ are the four complex U (A^)-valucd bifundamental scalars of ABJM; and A^, (// = 0, 1, 2) 
are worldvolume U{N) x U{N) gauge fields and e",/i" and are generators satisfying 



(E^)^ = £;-^, [£;-^,e"] =0, 



(4.2) 



Kab is the Cartan matrix and we have defined 

(Ai - iA2) 



A = 



Z = XI+IX2, 



D^Y^ = d^Y^ + iA^Y^ - iY'^A^ , 



(4.3) 



7V-1 



A=Y,Aah\ 

a=l 

4 



A=2 



G{z) is an arbitrary holomorphic function, C(, is a constant and xi,X2 are the two space worldvolume 
directions. The vortex solution then obtains by solving the above equations for the scalars and the 
gauge fields. There are two particularly simple cases: 



1. By taking the trivial value for the holomorphic function, namely G = 0, we can obtain an 
approximate solution. Define Ka = '^h^ab- Then yu = and the equation for becomes 
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55 In I yap = ^X^fe-^afebfeP) with A a constant. Consider moreover the case when Ka is inde- 
pendent of a, in which case we can choose a solution with g = ly^p independent of a, so that 
ddlng = AKag. Then near z = 0, the solution is approximately \y° 

JV-l 



,.a,2 _ gAK.I.p 



(4.4) 



a=l 



We will not calculate the gauge fields, since it is not clear how to interpret this approximate 
solution. 

2. The simplest exact solution is obtained for the simplest nontrivial G. In this case, taking 

,M 



y 



r 

Cn 



reduces the y"" equation to dd\n\y 



M|2 



and from the y equation we get 

/ N-l \ViV 



Vm 



= G{Z) n Ca 



(4.5) 



a=l 



Then, if we choose the simplest holomorphic function G{z) = z—zq, the equation ddln {yuf 
is identically satisfied and we have a complete solution. Specifically, 



Y^ 



>^ a=l 



+ e 




(4.6) 



= v\ ^4 = 2,3,4. 



For this solution can be interpreted as a supcrgraviton, we have to first argue that it corresponds 
to a classical pointlike object in spacetime. This is not too difficult to see. First notice that the 
corresponding coordinate Y^ is fixed and, according to the usual Matrix theory definitions, a VEV 
proportional to the identity corresponds to a fixed classical coordinate. Further, at the position of 
the vortex, z = zq, Y^ = Q is also fixed, so that this object is extended only in time and not in any 
of the parallel or transverse coordinates. 

Let's now calculate the gauge fields for this solution. To do so, we need to compute Aq,Aq,A 
and A (or equivalcntly, Ai, ^i; ^2)- Since F12 = F12, we can choose A = A. First, note that 
F12 = diA2 — d2Ai implies that dA = \[diAi + ^2^2 + ^-^12]- Consequently, in the Coulomb gauge, 
diAj, = diAi + ^2^2 = 0, the magnetic field takes the sussinct form B = F12 = AdA/i. Using the 
properties of the generators, 
we get 



Y'Y, 



Y^^=i |yaP^") so that substituting into (4.2) and integrating. 
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On the other hand, to solve for the we substitute the expression for the commutator [y^, y^^] into 

DoY^ = iAoY^ - iY^Ao = i=^v'^ Y\yI , (4.8) 

k I J 



to obtain 

AoY^ - y^io = ^^-"^(11 l^'^l' ) E r^l^ - -op/^ , (4.9) 



from which (after choosing Ao = 0) we can read off 

\b=l / a=l ' 

Note that for this vortex solution, the magnetic field B = F12 goes to zero at z = 2;o (the position of 
the vortex), but the electric field Fqz = —dAo diverges at that same location. 

This vortex solution looks like it could stand in for the spacetimc supcrgraviton, since it is a 
pointlike object in spacetime, carrying DO— brane charge. This last property is perhaps not obvious. 
However, as a massive, classical, pointlike spacetime object in ten dimensions this vortex must carry 
a charge corresponding to its momentum in the eleventh dimension and, as there is no other possible 
candidate, this must be identified with a DO— brane charge. As we don't have a dual description of 
this solution though, it is not easy to check this assertion explicitly. Wc will see that in the related 
case of the vortex of the massive deformation of ABJM that the corresponding object does indeed 
carry DO-brane charge. 

However, it is also easy to see that the solution has infinite energy, since the energy density of the 
magnetic field increases away from z = zq. Note also that the solution has the complex coordinate 
Y^ oc {z — zqY^^ which, in the large N limit exhibits a step function-like behaviour (0 at the vortex 
position and 1 away from it). In other words, away from the vortex, the solution represents a two- 
dimensional worldvolume growing at infinity, which could be identified with a 2-brane. As similar as 
they sound, this solution is different from Blonic branes of [27, 28], where a single coordinate X ~ 1/r 
signals a string extending to infinity at r = 0, as well as the self-intersecting M2— brane of [27, 29, 30], 
where the two complex coordinates^ s and t are related by s on c/t (c= constant). In each of those 
cases, a new brane or string "grows" at the position of the singularity, and represents the spacetime 
intersection of branes. Charge conservation then implies that the world vohimc flux must flow through 
a string or brane which must either extend to infinity, as for the Blon and self-intersecting M2, or 
end on yet another brane [27, 28]. This latter solution, ending at r = ro and finite X{ro), is half of 
the D ± Fl ± D brane configuration, and provides an example of charge conservation by ending on 
another brane with X{ro) finite and X'{ro) = 00. 

The vortex solution of pure ABJM is clearly an example of the latter. To be precise, since 
Y^{zo) = but y^ (zo) = 00, the only way to enforce spacetime charge conservation would be to 
have DO-branes at the endpoint. The situation is muddied however by the fact that Y^ grows in the 
direction of the original 2-brane (as z — ;> 00). 

In [31], it was demonstrated that the vortex solutions of pure ABJM described here have the same 
supcrsymmctrics as self- intersecting M2-branes. This fact alone likely means that wc can interpret 

^One of these is a worldvolume coordinate, t = xi + 1x2 while the other, s = X* + iX^, is transverse. 
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the solution away from zq as a self-intersecting M2-brane, with the caveat that DO-branes be added 
to the z = zo point to enforce charge conservation. It seems then that the vortex solution should be 
thought of as a bound state of self-interesecting 2-branes and DO-branes at zq. This interpretation 
matches nicely with the picture we will find in the massive ABJM case, where the vortex is a bound 
state of dielectric L>2-branes blown up into L>4-branes, and L>0-branes at the vortex position. The 
difference of course, is that in this case the 2-branes are infinite in extent, and have correspondingly 
infinite total energy. 

To summarize; even though the vortex solution looks like a supergraviton, it is hard to get actual 
physics from it. Essentially, this is because we will still need to have some appropriate regularization 
that to render the energy finite before we can subtract the 2-brane contribution and interpret the 
remainder as the supergraviton. There is no reason why this cannot be done in principle, we just 
have not found a simple way to do it. By contrast, it turns out to be much easier to work with the 
massive ABJM theory where everything is wonderfully finite. To this end then, in the next section 
we study the finite energy solutions of the massive ABJM theory. 

5 Massive ABJM solitons 

Following several earlier works on massive deformations of BLG theories [2, 32, 33], a maximally 
supersymmetric (A^ = 6) massive deformation of ABJM was recently proposed in [12]. However, it 
is not clear exactly what the brane interpretation of the deformed theory is. One suggestion is that 
before the brane backreaction, the background corresponds to the maximally supersymmetric type 
IIB pp-wave. In [26], a different conjecture, for the gravity dual of massive ABJM was proposed, 
based on a Zk orbifold of the LLM solution corresponding to the massive deformation of N M2-branes 
in flat space. 

For our purposes, we will need to understand the massive deformed background in which the N 
M2-branes move {before back-reaction) i.e. the analog of the C^/Zj^ orbifold for pure ABJM. Once 
we have computed the background dual to the deformed ABJM model, we proceed to analyze the 
various solutions of massive ABJM. From the point of view of the ABJM worldvolumc, these may 
be classified as brane-filling, vortex or kink type, according to their codimension (zero, two and one 
respectively). We will provide a brane interpretation of these solitons. 

5.1 Dual of massive ABJM deformation 

To begin, we will compute the spacetime background corresponding to the massive deformation of 
ABJM. Recall that in the case of pure ABJM, the spacetime background is set up by N M2-branes 
living at the tip of the cone on i?^'-*^ x C^/Zj;. We will deal with the Zj^ orbifolding at the end, but for 
the moment we would like to understand the background that replaces the flat eleven-dimensional 
i?^'^ X in the presence of the massive deformation. To this end, following the suggestion in [12] 
that it could be related to the IIB maximally supersymmetric pp wave, we consider the k = 1 case 
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first and notice that the type IIB gravitational wave 

ds^ = -dt^ + dx^ + {H - l){dx + dtf + dx^ , (•5 ]^) 

e't' = gs, B = 0, 

is T-duaF to the fundamental string (Fl) solution 

ds%^ = H-^{-dt^ + dx^) + dx^ , 
Bt^ = l-H-\ (5.2) 

where H = 1 + Q/r^ is a harmonic function of the transverse coordinates x and the Fl solution is 
such that all fields are trivial at x — >■ 00. In turn, using 

dslf = e^'^/^ (dslo + e^'l'{A^dx'' + dxuf^ , 

F^tll = e^^/=^(F„6c<i + 4^[ai^M) ' (5-3) 

^abcll ~ ^o-bc ) 

the -Fl solution lifts to the membrane (M2) solution 

dsl = H-^/^df + H^/^df, 
F4 = dy^ A dy^ A dy"^ A dH-^ , (5.4) 
Q 

in M-theory. Here we define {y^,y^,y'^) = {t,x,y = xu). On the other hand, the maximally 
supersymmetric type IIB pp-wave , 

ds^ = 2dx^ dx~ — iJ?x^ [dx'^^'^ + dx^ , 

-^+1234 = -^^+5678 = M ) (5-5) 

= g,; B = 0, 

can be understood as the gravitational wave obtained in the presence of the constant flux -F+1234 = 
^+5678 = A*) which also modifies H—1 from Q/r^ to — /x^r^. Correspondingly, the IIA solution T-dual 
to this pp-wave is found by using the full IIB — )• IIA Buscher rules (including those for the RR sector) : 

io« = |B<f + 2^, (5.6) 

A, = -5jf+aBo,; io = 0, 



'^The easiest way to see this is to use Buscher's T-duality rules [34, 35, 36] which, in the absence of RR fields are the 
same for IIA ->■ IIB and IIB — >■ IIA, namely: 

1 ~ Boi _ goigoj — BoiBoj 

900 = — ; goi = — ; gij = gij 

500 300 900 

i>()i ~ , iiij = tiij H 

.900 

4> = (p- - In 500 
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to get 



str 



H-^ {-dt^ + dx^) + dx 




(5.7) 



-^1234 



-^5678 — M ) 



H 



1 2 ->2 
1 — H X . 



This can be interpreted as a fundamental string in the background of the constant flux -F1234 = 
-^^5678 = which, again, modifies the harmonic function. Finally, this lifts to the M-theory solution 



F4 = dy^ A dy^ A dy^ A dH~^ + ii{dx^ A dx^ A dx^ A dx'^ + dx^ A dx^ A dx'^ A dx^) , (5.8) 



This can again be interpreted as an M2-brane in the presence of constant transverse flux F1234 = 
-^5678 = M- This flux (in particular |-Fp) means that H now satisfies d^H = so that, again, 

the harmonic function is modified. There are two points about this geometry that should be noted: 
before the T-dualizing, we could always rescale ji away by sending x+ — >■ x+//x and x" jJLX" . After 
the T-duality and M-theory lift, this is no longer possible. Now the scale jj, has physical meaning. 
Secondly, while there is nothing particularly interesting about the point \x\ = in the type IIB 
metric, here it is potentially singular. 

This solution obtained can be interpreted as an M2-brane with a constant flux preserving the 
same manifest i?-symmetry as the massive deformation of the ABJM field theory. Specifically, the 
SU{4:)r is broken to SU{2) x SU{2), by the splitting of the scalars into the 1,2,3&4 and 5,6,7&8 
directions. This gives us confidence we are on the right track. For this, and other reasons given in 
[12], we can say with some confidence that we have found the spacetime background corresponding to 
the massive deformation, at least in the k = 1 case. For general k, we still have to apply the action 
(inherited from the pure ABJM case) on the transverse coordinates. Since it acts by — >■ e^'^*/*^y^ 
the solution remains intact and the background is valid at any k. 

There is one slight subtlety related to the amount of symmetry preserved under T-duality in 
the above construction. We started with maximal supersymmetry (32 supercharges, or AA=16 in 
three-dimensions), and ended up with a solution that cannot have as much since, in M-theory, the 
unique backgrounds with maximal susy are flat space, AdS^ x S^, AdS-j x 84^ and the maximally 
supersymmetric eleven-dimensional pp-wave that is their Penrose limit. Since our solution is of the 
M2-brane-type, we must have that at least the constraint T^T^V'^il) = tp hold. This reduces 32 
supercharges to 16 (or = 8 in three-dimensions). In addition, the presence of the transverse 
flux -F1234 = i^5678 = M further reduces the number of supercharges to 12 (giving AA = 6 in three- 
dimensions). Naively, this breaks the i?-symmetry to an SU{2) x SU{2) x U{1) which rotates the 
transverse coordinates Y^. The problem is that in [12] it was argued that the massive deformation of 
the ABJM model retains the full SO{6) i?-symmetry of N=6 in three dimensions, so the same should 
apply here. 

On the other hand, in [32], it was argued also that in the BLG (or, equivalently, the A=2 ABJM) 
case, the massive deformation with M2 on Rf xT^ is equivalent to the IIB pp-wave. There, the massive 
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BLG model was observed to satisfy 32 super symmetries, only 16 of which are linearly realized, and the 
remaining 16 are nonlinear. The fuzzy sphere background breaks the 16 nonlinear supersymmetries, 
giving 16 genuine supersymmetries. These fuzzy sphere vacua correspond to D3-brane giant gravitons 
in the IIB pp-wave background, which preserve the same 16 supersymmetries in addition to satisfying 
the same zero energy condition. 

Returning to the case at hand, it seems that the solution to the problem of supcrsymmetry under 
T-duality is that one has to consider the full string theory (or rather Matrix theory) and not just 
the supergravity background. Then one needs to consider DO-branes in IIB in the presence of the 
pp-wave. This preserves a total of 16 supercharges, instead of the full 32. After T-dualizing, in 
the M-theory we have an M2-brane with constant transverse flux and H = 1 — n^x^, on top of 
the ABJM Af2-brancs. This configuration again preserves 16 of the supercharges, so there is no 
problem. Of course, the above argument applies just for the k=l &; 2 ABJM model (that has no 3/4 
supersymmetry reduction), otherwise we need to take a quotient of the target space. 



5.2 Brane-filling solution: the fuzzy sphere 

The first type of solutions we will analyze in some detail are brane-filling, i.e. ground states, which 
can be interpreted as branes with a transverse extension. In [12] it was shown that the maximally 
supersymmetric ground state of the massive deformation of ABJM is given by an equation associated 
with a fuzzy sphere. In particular, the four complex scalars of ABJM as split as = (Q", i?") and 

the solution is given by = and = = y^, with 

-G" = G^gIg" - G^G^G^ . (5.9) 
The irreducible matrices that solve eq.(5.9) - that we will call the GRVV algebra - are 

{G^)m,n = \/rn--T Sm,n , 

{&)m,n = \/{N -m) 6ra+\,n , (5.10) 

iGi)m,n = Vm - 1 6m,n , 

{Gl)m,n = \/iN -n) Sn+l,m ■ 

Based partly on the fact that G'^Gt = iV — 1 for this irrep, this solution was conjectured to represent 
a fuzzy S^, but this was shown to not be the case in [6]. It is instead a fuzzy S^, as could be guessed 
by the fact that G^ = G\ (so that the scalar represented by the imaginary part of G^ is fixed to zero). 
Moreover, in [7] it was shown that if we mod out by the U {N) gauge transformations that leave the 
GRVV algebra invariant, the resulting bifundamental matrices can be viewed as fuzzy versions of 
Killing spinors on the S"^. The action for fluctuations around the solution is then the supersymmetric 
Z)4-brane action, compactified on S^. The fields can be expanded in fuzzy spherical harmonics made 
up of either Jj or Jj, 

J, = {ar(,G^Gi = {a,r^jP^^{a,)^"jf'a, 

(5.11) 

Ji = {^iY rG^^G^ = {^iY aJa^ = {(^i)p°'Ja^ , 
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in the same way as the regular spherical harmonics are built from Euclidean coordinates Xi. In the 
classical limit, both the Jj and Jj tend to the same Xi, up to a normalization constant. The D4- 
brane compactified on the S'^ ground state then corresponds to a DS-hiane giant graviton wrapping 
a 3-sphere 

\Z^f + \Z^f = , (5.12) 

in the maximally supersymmetric type IIB pp-wave. Here the 4+4 coordinates transverse to the IIB 
pp-wave are = (Z", Z") with 

Z^ = X^ + iX'^ , Z^ = X^ + iX^ , 

Z^ = X^ + iX^ , Z^ =X'^ + iX^ . 

We also have to mod out by Zk which acts on the transverse coordinates as Z^ e^'^^^'^Z^, and Z^ 
corresponds to in M-theory. 

5.3 Vortex solutions 

Moving up in co-dimension, in [24], a BPS vortex solution was found for another massive deformation 
of ABJM. This time the deformation is given by the superpotential term SW = fj,Tr[2'^yVa], and 
admits the supersymmetric ground state 

= VF'f" = (5.14) 
V 47r 

Here we have written = (Z",Wt") instead of {Q°',R'^), since the different mass deformation 
implies a different kind of split. While this deformation manifestly preserves an SU(2)r symmetry 
and N = 2 supersymmetry, it was demonstrated in [25] that it in fact preserves the full A/" = 6, and 
it matches the GRVV deformation upon a redefinition of the fields. The BPS vortex solution is given 
by the ansatz 

Z" = W^"" = f{x)G'' , 
= a^(x)G"GL 

= a^(x)GtG", (5.15) 
f{r,e) = ^g{r)e'-\ 
ai{r) = eij^[a{r) -n], i = l,2, 
where o(r) and g{r) are real functions satisfying 



a(r) = =frg'{r)/g{r) , 
{\ng')" + l{lng-'y+4n'g\l-g') = 0, 



(5.16) 



together with the boundary conditions a(0) — n = 0, ng{0) = and g{oo) = 1, a(oo) = 0. The 
vortex has magnetic flux $ = —kiJ,n/2 and zero angular momentum and, at infinity it goes over to 
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the GRVV ground state, i.e. the fuzzy S'^. Consequently, at the position of the vortex r = 0, and 
Z°' = H^t" = 0. It has no extension in the transverse coordinates either (the fuzzy sphere shrinks to 
zero there), and so can be interpreted as a point in both worldvolume and transverse space directions. 
This implies it can be identified with an object carrying DO-brane charge^ [26]. The U(l) symmetry 
under which the vortex is charged was identified in the dual gravity background - the reduction of 
the LLM solution for N M2's in flat space - as J = kQo + NQ4 where Qq and Q4 are the DO-brane 
and D4-brane charges respectively. One interpretation of this solution is as a bound state of an object 
with DO-brane charge k, the vortex, with the D4-brane wrapped on the fuzzy 5^ corresponding to 
the ground state. This picture was confirmed through a DO-brane probe analysis with the mass of 
the probe at the minimum matching the soliton mass of k^. 

Evidently, the vortex truly docs represent a charge k DO-brane object which, in the large k classical 
limit becomes a classical spacetime object, a supcrgraviton. It also follows that the vortex solution 
of the massive ABJM theory is indeed a crucial ingredient in a search for a Matrix model of ABJM. 

We note that if the proposal of this paper is correct, we should also find ABJM solutions corre- 
sponding to branes, at least D2 and D4, that are completely transverse, i.e. that have a worldvolume 
in the transverse scalar directions, but are pointlike on the worldvolume. These would be vortex 
solutions with fuzzy 2- and 4-dimensional spaces respectively, at the position of the vortices. We 
have, unfortunately, not succeeded in identifying such solutions. 

5.4 Kink solutions 

According to our proposed Matrix theory-type interpretation for ABJM, we also expect to find ABJM 

solutions that can be interpreted as 2-brancs with one spatial worldvolume coordinate parallel to the 
M2-brane and one perpendicular to it. Such a solution would be a kink soliton of ABJM model, which 
at the position of the kink can be interpreted as a "fuzzy circle". The problem however is two-fold: 
first, we need to find a kink solution of the ABJM model and second, we have to understand what 
exactly such a "fuzzy circle" is. Odd-dimensional fuzzy n-spheres have proven notoriously subtle to 
construct in the past when n>3 [37]. The fuzzy circle, on the other hand, appears not to even have 
been constructed yet! We will see however that in a particular limit, the fuzzy sphere plays the role 
of a fuzzy circle.^ 

The general problem of finding a fuzzy circle can be readily described. For a fuzzy sphere, the fuzzy 
spherical harmonics constitute a complete subset of the set oi N xN matrices, with a maximal angular 
momentum of N'^ — 1. This also means that the composition of spherical harmonics, Y^Yi ^ ^ki^n 
is still valid in the fuzzy case as long as A; -|- Z < L^ax = — 1- For a fuzzy circle, the "fuzzy 
spherical harmonics" (or "fuzzy Fourier modes" ) would also need to satisfy similar completeness and 
composition relations. These, however, appear to be quite difficult to obtain in any other way than 
as a limiting case of the fuzzy sphere. For instance, completeness requires that the elements of 

*In [25] and [26], more general vortex solutions were given, based on a general discussion of the BPS equations of the 
GRVV massive deformation. Because of the identification of the massive deformations up to field redefinitions however, 
and the fact that vortices carry topological charge, all the vortex solutions must represent the same object. This soliton 
has charge k and mass k/i. 

different limit was considered in [38], where a fuzzy circle arose from a "fuzzy cylinder" construction. We thank 
Yolanda Lozano for pointing this out to us, after the first version of our paper was posted on the electronic archives. 
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the N X N matrices need to be arranged as circle Fourier modes of increasing order n which in turn 
makes composition hard to accomphsh. 

We now turn to studying kink solutions of the massive AB JM Lagrangian. Naively, it might seem 
like a good place to start is with the construction of [24, 39] of the half supersymmetric kink 



r(^i) 



k fi 



(5.17) 



47r 1 + e-2M2:i 



However, at the location of the kink, xi = and Z" = W^°' = /xG", resulting in a degenerate fuzzy 
S'^ X S'^. This is not the kind of solution wc arc looking for^'^. 

In search of the fuzzy circle kink solution, we will look for solutions with no gauge fields, and half 
the scalars set to zero. In the notation of GRVV [12], this corresponds to exciting only the of 
= The action on such a solution would be 



S= d'x 



+ 



fiQ' + ^{Q'QlQ^-Q^QW) 



(5.18) 



We can readily check that the ansatz = f{xi)G^, 



-fik 
"2^ 



G^, where G\G^ are matrices 



satisfying the GRVV fuzzy sphere algebra, results in the following quartic action for the function /, 

2" 



S=^ I dxi 



(5.19) 



The resulting equations of motion are solved by the usual kink solution f{xi) = tanhxi. However, 
we note that only half of the equations of motion of the full action eq.(5.18) are solved by this ansatz. 
This "near miss" suggests another avenue: Completing the square in the action eq.(5.18) produces 
the BPS equations 



dsQ" 



k 

27r 



(5.20) 



MQ' + ^(QVig'-QViQ'), 



whose solutions will be the half supersymmetric solitons of the theory. The ansatz = /T^; 
gT^, where the are matrices that satisfy the algebra 



fJ-^2 fJ-^f rj-il rp\ rj-\\ rJ-^1 



(5.21) 



^"One could perhaps interpret this solution as a 5-brane with 4 directions transverse to the ABJM worldvolume, but 
this is less clear 
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collapses the BPS equations to two coupled differential equations for / and g whose solution is 
/(^i) = 



-Ci 



1 - 27rQ;e-2CiC2/kexp(^e2Mxi) 



g-ClC2/kg^p 



, 2^X1 

2/ik 



(5.22) 



27r y 1 - 27rQe-2CiC2/k exp(^e2/^^i ) 



After fixing one of the integration constants as Ci = — fc^/27r we get 

k 1 



9 



f 



27r 



l_4^^2exp(-^e2M-) 

.4 



(5.23) 



y/l-^A2exp(-^e2M-) 



exp 



47r/i 



This solution has the following properties: As xi — >■ ±00, / — > 0. Moreover, /(xi) has a maximum at 



27r/x 



47r^ 
fee 



of value 



frr 



A 



L 



ke 



+ 1 



a 



(- 



ke 



a 



(5.24) 



(5.25) 



(5.26) 



where -Z^(^) is a Lambert W function that satisfies the Lambert equation, 

,_, dL ^ ,1 
z{l + L)—=L, z + — . 
az e 

This solution can be interpreted as a (non-topological) kink with regard to the field /, while nothing 
particularly interesting happens for the field g. As for the transverse fuzzy circle, note that the 
can be taken to be 

t1=G^ ^2 = V^G^ (5.27) 

and hence define deformed fuzzy spheres. Up to now, the introduction of the parameter a. seems to 
be superfluous. However, note that if we take the a — )• 0, then what we have is a BPS solution with 
r-*^ = and = 0, yet satisfying the nontrivial relations"*^^ 



rplrp\ rpX rp\_rj-}\ rp'2 



0, 



(5.28) 



This is ample justification for our associating them with the fuzzy spherical harmonics of the fuzzy 
circle. In this limit, the solution becomes 

k 



9 



27r 



ylexp 



(5.29) 



lixi 



-,2/ixi 



47r/i 



^^Of course, if = from the start, these equations are trivial. They are only nontrivial on taking the limit of 
Q ^ 0. 
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and is still a kink, since / has the same asymptotics with a maximum now at 



while the field g exhibits a simple monotonic increase. The interpretation of this result is that we have 
taken a fuzzy sphere and shrunk one of its directions (G^) to zero, thus obtaining a fuzzy circle-*^^. 



6 Graviton scattering in massive ABJM background 

There are two cases of interest for the scattering of gravitons in the background given in eqs.(5.8). 
These correspond to the two types of vortex scattering in the massive ABJM model that we want to 
match to. There, the two relevant cases are when the vortices are separated in the directions parallel 
to the ABJM worldvolume, or in the transverse (scalar) directions. 

In one case the gravitons are separated only in the two directions parallel to the ABJM world- 
volume, and the 11th dimension is transverse to the worldvolume. In the second, they are separated 
only in the directions transverse to the ABJM worldvolume, and the 11th dimension is parallel to it. 
This latter case is a more direct analog of the BFSS analysis which we reviewed previously, whereas 
the former case is new. As in the BFSS model, we will employ 't Hooft's trick of replacing the 
graviton-graviton scattering with the scattering of a graviton Shockwave with a probe wave, but here, 
it is not clear a priori that it will give the complete result. 

As in the BFSS case, the graviton Shockwaves satisfy the wavefunction equation, 

Aig,h..=QS^{x,y). (6.1) 

We will content ourselves with some cursory comments about the solutions of this equation here and 
leave a detailed analysis of the solutions for the Appendix. As earlier, we will denote the worldvolume 
coordinates by y (with |y| = y), and transverse coordinates by x (with, correspondingly, |x| = r). In 
the transverse separation case, we find 

h-^iy = 0,r)^C'^, (6.2) 

as in the BFSS case, but since the result is independent of /x, we could have imposed the wrong 
boundary condition, so we will refrain from using this result in the following. In the case of parallel 
separation on the other hand, we find that 

h—{y,r = 0)^Ce-"'^\y\. (6.3) 

The next step in this computation is to scatter a scalar probe-graviton off the source-graviton. Since 
it is clearer how to deal with the transverse case, we will do this first, more as a guide than anything 
else. 

^^As a final point, note that our choice of representation for the G's is not unique. We could take for and either 
the irreducible GRVV matrices (as we have), or the reducible matrices = Sm+N/2,n, G\ = 5m,n+N/2 and G^ = Smn 
with m,n< N/2. 
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1. Transverse separation: Expanding the D'Alembertian in the probe graviton equation 

= 0, (6.4) 



□4f + d+h++d. 



gives 



9^ 



-{do - mf 



(6.5) 



As in the, perhaps more famihar, BFSS case the momentum-space interaction potential is computed 
through the S-matrix 

(6.6) 



S'''"\pin - Pout) = 

where 0^°^ solves the unperturbed equation 



</>iS(.)*F(z)<A;°)(z), 



(6.7) 



If we factorize the plane wave 0(0) as (/.W = ^"^;^"^V (x), define 



Po 



and E = (po)'^ /{^/ipo), we find that the function / satisfies 



dl-f + 2E / = 0, 



(6.8) 



which is nothing but the equation of a d-dimensional harmonic oscillator. Therefore, if we were 
interested in solutions that drop off at infinity, the solution would be given in terms of Hermite 
polynomials as 



(6.9) 



E 



with the quantization condition E-i = nj + 1/2. 

However, we have \x\ < and, as explained in the Appendix, we need to impose Neumann 
boundary conditions there instead. Consequently, the "plane wave solution" of (6.8) is 



/ 

fpo,i 
2npo 



= e 



El ^ 

2 ' 2 



i 



{P0,i? 
2fipo 



(6.10) 



Here the relative normalization C1/C2 could be defined by the Neumann boundary condition at 
r = 1//X, but that would result in two constants for each dimension. We will see in what follows that 
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another relative normalization is preferred. Note that because = X^j(po,i)^' Po can be called the 
(plane wave) momentum, and we will fix the normalization by requiring that we obtain the fiat space 
plane waves in some limit. 

Indeed, from Appendix C of [40], we know that the harmonic oscillator wavefunctions become 
plane waves as, for example n — >■ oo or mu — >■ 0. Explicitly, 



<t)n{x) 



^!l^:i)!l_g-m.xV2^ ( V^x) ^-^i ^1 cos (2^/^x) (6.11) 
(0F2"n!)^/^ ' ^/^ \ n J 



where kn = y/2mu)n = I'KnjL. This limit, together with the relations 



(6.12) 



(2n + 1)! ^ ^ f 3 2 



H2n+i{x) = (-1)"2 ^, a; iFi ( -n,-;x 



generalize to the limit of our solutions, giving plane waves for large po^i/fi. We fix the relative 
constants Cii/C2i such that we are able to sum cos(- • • ) + «sin(- • • ) = e*^"') after the limit. As we 
will argue further on, knowledge of the precise values of these constants is not necessary, only that 
the plane wave limit exists. The overall constant in front of the solution follows from a relativistic 
normalization. 

We now reverse the initial logic, take the po,i to be independent, and instead constrain E by the 
relation + m? — Py = = J2i{Po,i)^- The interaction potential is then obtained from 



S'^^\E; Pin - Pout) = J dxodyd 



^out:py,m,pQ 



/i++(y,x) 



(9o ~ fn 



tn:py,m,po 



A\fW^^E°^ 

= 27r5{E°-'-E^^)V^i{E;pin-Paut). (6.13) 

with the understanding that we need to make a Fourier transform back to x space at the end. However 
note that since /po(x) ^ e*P°^, it is not even guaranteed that we get a function of \Spo\ = — pg"*! 
after the x integration. Assuming that it is nevertheless such a function, after the Fourier transform 
back to position space we find that 



^ a}^! d^po e-'^'P" J d«f :^tU^/i++(y,r)/|„(x) (6.14) 
In particular, when y = 0, 

= 0,z^ = ^ 1 d'Spo e--^^-° I d'xJ^Mh-^^iy = 0,r)/|„(x) (6.15) 

While it is not clear whether the integral is dominated by the low r region, what is clear is that we 
cannot have a result depending on the dpQ unless we are in the plane wave limit. This is obtained in 
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the limit of high momentum poi- seems then that only if we probe with sufficiently high momentum 
we will have a well-defined interaction potential. 

In any case, the outcome is that we need to replace the functions fp^ with plane waves and obtain, 
as in the BFSS flat space case, 

,,4 



V^l{y = 0,z) 



mv 



= 0,r). 



(6.16) 



2. Parallel separation: 

To write down the equation for a probe graviton, we need to first define the coordinates and metric 
more a little more carefully. We start with the metric 



ds'' = H-'''\x){-de + df) + Hy^dx\ 
dx^ = dr'^ + r'^dn'^; x = {x^ , ...,x^) , 



.17) 



and dimensionally reduce on xn. Operationally, we choose xn as the fibre direction of the Hopf 
fibration of S'^ over CP^ which acts on the Euclidean coordinates, Y"^, on the 5^ through the overall 
phase = e^^iiy^. Since the coordinates satisfy |y^p/r^ = 1, this means that if we write 
X = {x',xii), then r = \x\ = \x'\ and dx^ = \dY^\^ = |dy"^p + dx^^ becomes dx'^ + dx^i- Since the 
metric is translationally invariant in Xn, we can schematically rewrite 



ds'^ ~ i?-2/3(f ) (-dt^ + H{x')dxl^ + df) + H^/^{x')d^ 



(6.18) 



With this metric, the equation for the probe graviton takes a form similar to the transverse separation 
case, 



d' 



-{do - m) 



where A (x',m^) is obtained from 9| by reducing on xn, and setting ^(xn,.. 



(6.19) 
s^P^ii (/)(...) and 

p = im. Noting that we could rewrite —dt'^ + H {x')dxii = 2dx^dx^ only on slices of (nearly) constant 
x' , we picked the simplest choice x' ~ and fixed it so that we can form the x"*" combination for 
defining without problem. 

Now we can proceed as with the transverse case. Specifically, the momentum-space interaction 
potential obtains from 

S^'Hpin - Pout) = J d''z <AiS(^)*y(^)0i?(z) , (6.20) 

where (f)^^^ satisfies the unperturbed equation associated to eq.(6.19). Here again, a plane wave ansatz 
of the form = ^^^^^==— ^/(x) will solve the equation provided 



We again obtain 
S'^'^\E;pin-pout) 



I 



[H{x) -fy)+A (f, m^)] fix) = 

/i++(y,x 



(6.21) 



dxo d y d X 



out:py,m,po 



(do - m) 



m:py,m,po 



(6.22) 



5[E'^^ - E^'^){E''^ - mf 



J d'y e^y(^T-PT*) I d'xf^ih++iy,r)f^JL 
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or, after manipulations similar to the transverse separation case, 



{E 



m 



4E 



d'Spo I dJxf;p{xrh-^+{y,r)qniS) 



(6.23) 



We have, of course, already noted that equations (6.19) through (6.21) only make sense if x 2± is 
implied, in order to be able to form the combination. This could be alleviated by integrating x 
only in a neighbourhood of zero, in which case 

{E - mf 



AE 



-h-^^{y,r = 0) 



/ d'^dpo [ 

J J\x\c^O 



d X f ^y^t(x) f pin{x) 



(6.24) 



In both cases however, we still have the same problem that we observed in the transverse separation 
case. Namely, in order for the momentum space result to only be a function of 5po = p^^ — Pq^'', we 
need to only consider high momentum wavefunctions, which are effectively just plane waves. In that 
case, the integrations disappear, irrespective of whether or not x ~ 0, and we obtain 

XE-m'^'^ 



4E 



-/i++(y,r = 0). 



or, since h'^'^ oc e 



(6.25) 
(6.26) 



So much for the graviton-graviton scattering. 



7 Vortex scattering in massive ABJM 

In this section-*^^ we analyze vortex scattering in the N = 2 {ot equivalently, the U{2) xU(2)) massive 
ABJM model. The general N case is quite complicated, so we will not attempt it here. We will employ 
the vortex solution in the form of [25], as it is both guaranteed to be the most general solution^^, 
and also because the specific form of the solution will allow the use of previous known results for the 
Abelian-Higgs model and its Nielsen-Olesen vortex solution. The vacuum for the massive U{2)xU (2) 
ABJM model is 

yi^./^G-i; Y^ = \f^G''; Y^ = Y^ = 0, (7.1) 
V vr V vr 

so that, keeping = = and using complex notation for z = x^ + ix"^ means that in terms of 

= ay" + zA^y" - lY'^A^; p^y" = ay" + i^^y" - iy"i^ , 



(7.2) 



A, = ]^{Ai-iA2)] A-, = ]^{Ai+iA2); = ^(ii - zia); = ^(ii + zi2) , 



the BPS equations read 

Pjy^ = 0; ViY"^ = V2Y^ = , 

VoY^ - + f,Y^) = 0; VoY^+i{Pl^ + iJ,Y^) = 0, (7.3) 

AeM'^pp = i(y"(D'^y")t - D'^y^yJ;); Ae^^^p^ = z((D'^y")ty° - Y^D^'Y'') 

2iT 27r 



^^This section was done in collaboration with Toshiaki Fujimori. We thank him for allowing its use in this paper, 
^''since it was not found using any specific ansatz, but rather by analyzing the energy functional on a case-by-case 
basis. 
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where a = 1,2 and 



4tt 



T k 

The multi-vortex solution can then be written as 



(7.4) 



Y 



An 



kfi (O 1 



1 ddip 







27r \ 

io = 



^2 ^ /fc^ / 

V 27r I 1 



(7.5) 



2-/0 



Az — Az 





pV'. 



Here i?o(-2^) is an arbitrary polynomial 



Ho{z) = l[iz-Zi), 



i=l 



and the function real 'ip{z, z) is determined through the equation 



(7.6) 



(7.7) 



with the boundary condition at \z\ — t- cxd requiring Tp — t- log|i7o(^)P, and where Zi are the position 
moduli for n vortices. Correspondingly, the energy of this solution is nkji. Note that this is the same 
equation governing the vortices in the Abelian-Higgs model, so we expect that the same effective 
action governs their scattering as well. This is indeed the case. 

To obtain the effective action for vortex scattering, we let the parameters become functions of 
t, so that the fields of the static solution become a function of {z,z, Zi(t),Zi{t)). Of course, if we do 
that, the equations of motion are no longer satisfied, but we can still solve them order-by-order in 
time derivatives df ^ Zi{t). We write 



Y^{z, z, t) = Y(^) {z, z, Zi{t),Zi{t)) + y^?^ + ¥(% + .. 



(1) ^ ^ (2) 



(7.8) 



and similarly for A^,A^. Because they are multiplied by the zero-th order equations of motion, the 
second order fields do not contribute to the effective action, so we can stop at the first order solution. 
At that level (on-shell), 

^e// = j d^xC{Y^,A^,A^). (7.9) 
The first order solution is found to be 



y-a _ n. 4(1) _ 4(1) 
-1(1) — u, ^0 ~ ^0 










-,{z'di - z'-di)'ip 



(7.10) 



41) ^ 



^ z'didil^ 



A 



(1) _ 














Finally, substituting this solution into the Lagrangian, produces an effective Lagrangian for the moduli 



L 



eff 



1 



1 



ddil) H — z^z^ I didjil) H — ^(didi/jdjdijj — ddipdidjip) 



(7.11) 
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where we have used the equation of motion (7.7) and its boundary condition at infinity. This effective 
Lagrangian is exactly the same as the one obtained in the Abehan Higgs model for the Nielsen-Olesen 

vortex! So both the equation of motion (7.7) and the effective action for the nonabelian Chern-Simons 
ABJM model give the same result as the ones for the Abelian Higgs model. We can then rewrite the 
effective Lagrangian for large vortex separation as in the Abelian Higgs model as 

Le//~-A;/in + ^^|if -A:/xg^i^o(2/x|/-z^|)|i^-i^f , (g ~ 1.71) (7.12) 

1=1 i>j 

which means that at large separation, the interaction potential becomes 

Vint - -kfiqKoi2fiyy ~ -k^qy-——v'^ (7.13) 



8 Towards a Matrix model of massive ABJM? 



Since we have now computed the scattering of both supergravitons and vortices in ABJM, we can 
compare them to see if they match. Assuming that they play the role of supergravitons^^, the effective 
potential for vortices was found to be oc v^e~'^^y , with y the parallel separation. The v'^ dependence 
is a general characteristic of lowest order soliton scattering, and the e"^''^ = e~"^°'^ arises from the 
presence of a lowest-mass excitation in the "broken" phase, i.e. the fuzzy sphere background at 
infinity. By contrast, we have found that the supergraviton result is proportional to v^e""*!^. Here, 
the v'^ was due to the same calculation as in the BFSS case while mi ~ 9.1// appeared because of the 
necessary Neumann condition at r = l/n which, in effect compactifies the transverse direction, with 
mi being the lowest KK-mode. Unfortunately then, the interaction potentials for separation in the 
directions parallel to the ABJM worldvolume do not seem to match. 

One possible reason for this mismatch - assuming, of course, that our identification of the vortices 
as supergravitons is at least correct - appears when we realize that the same dependence appeared 
on the gravity side as in the BFSS model. Unlike in BFSS where the first interaction term is a 
one-loop effect, however, here on the matrix side we computed a classical effect. Clearly then, we 
need to think a little harder about what terms are supposed to match on both sides, and try to find 
a corresponding gravity term for the vortex calculation we did, and a corresponding vortex term for 
the gravity calculation. 

There is potentially one other subtlety, though a priori, it seems to give only small corrections, 
namely that the gravity background we took is not complete. The pure ABJM field theory corresponds 
to the IR of M2-branes on C^/Z^, but we also need to consider an appropriate vacuum. For pure 
ABJM, the vacua are just VEV's which correspond to a given position in C^/Z^. On the other 
hand, in the massive ABJM theory the same VEVs, corresponding to positions in the dual gravity 
background eqs.(5.8), are no longer solutions. Instead, the vacua are fuzzy spheres of radius r = l//i, 
giving an M5-brane wrapped on S^/Z^, which will backreact and modify the background, maybe 
smoothing out the r = 1//U potential singularity. 

^^The background of D4 on fuzzy is common to both vortices in the scattering problem and thus does not contribute 
to the interaction 
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Before we continue trying to fix the mismatch, it will be useful to build some intuition for the 
problem by looking in more detail at the matching in the BFSS model. 



8.1 BFSS matching: details 

The review of the systematics of Matrix theory matching in this section will follow closely [41]. With 
the rescalings r = u/R and X* = y^/Mp, the bosonic part of the BFSS Matrix theory action (2.1) 
becomes 

making explicit the fact that the coupling of the theory is Mp (or that the loop-counting parameter 
is Mp). At L loops then, the BFSS effective action takes the schematic form 

Sl = Mf-"" J du h {y\ D^) = RMf-'' J drfL (M|X\ R-^D-,) , (8.2) 

and, as advertised in section 3, there are two sources of Mp-dependence in the interaction potential; 
one coming from the coupling constant and the other from the use of Mp to build a unit of length. 
Further analysis, including keeping only Lorentz invariant terms, gives the result 

S. = yd.W|(M,.)-"(^^,jj^^) (8.3) 

SO that, for the first few terms 

4 6 
Cii V Ci2 V 

'^MfR^V^ ^ M}i>R^ ^ ^ ■■■ (g_4) 

4 6 
C21 V C22 V 

+Ml2i?3 + ^^8^5 ^04 + - 
r = n I ^31 C32 

^ MfR^ ^ MfR^ ^17 ^ - 

Notice that the diagonal terms in this series that contribute to the classical general relativity scattering 
themselves receive contributions from all Matrix theory loops. The Einstein action has an infinite 
number of vertices, thus tree level scattering gives an infinite series, but they all come with integer 
powers of Mp ~ 1/Gii ~ 1/(kii)^. On the other hand, the off-diagonal terms in the Matrix theory 
loop expansion must come from quantum correction to the effective action of gravity, since they have 
nonintcgcr powers of k^j^. 

As an example, a source-probe approximation for the true classical general relativity interaction 
potential of gravitons gives 

N2 2 15 A^iA^2 225 NIN2 „ / 



^ 2r'' ^ 16 Rml r7 + 64 R^Mf r^^ ^ ^ Vr'V ' ^^'^"^ 

where the result is linear in the probe momentum N2/R because of the approximation, and Ni/R is 
the source momentum. The symmetrization of the last term matches against a 2-loop Matrix theory 



26 



result 

225 1 NiNl + N^N2 

32 Rm},^ ri4 2 ■ ^ ' 

In any case, we see that in BFSS, all the dependence on v, r, Mp and R could be determined from 
general principles. It seems pertinent to ask if the same could not be done for ABJM? The answer it 
seems is no, at least not for the pure ABJM theory which is conformal. There, all Afp-dependence 
is tied to the r-dependence so we cannot take advantage of having a dimensionful rescaled field and 
coupling, like and Mp, to use in dimensional analysis. Perhaps in the massive ABJM that we 
used for calculations something can be done since, for instance, the size of the fuzzy sphere vacuum 
/ = y^kfj,/{4:TT) sets a dimensionful scale. This remains an open issue. Sadly then, it appears that 
the intuition of BFSS matching is not particularly applicable to our case. 



8.2 ABJM matching and effective Lagrangian 

Let's step back for a minute and take stock. Our analysis up to this point has not succeeded in 
producing a concrete Matrix theory for the ABJM model. However, we have learnt a significant 
amount about the structure of the problem that we are in a position to speculate on what the sought 
after Matrix theory and its matching to gravity would look like. To outline the way forward as clearly 
as possible, we will focus our speculations on three specific points; the matching of parameters on 
both sides, what calculations can be done on the ABJM gauge theory, and what on the gravity side. 



• Parameters: The parameters available on the gravity side of the calculation (for the back- 
ground dual to the massive deformation of ABJM) are: the mass deformation fi, the inverse size 
of the 11th direction 1/R, the Planck scale Mp, the source and probe 11th momenta Trig = Ng/R 

and rup = Np/R and mi, the effective mass of the graviton wavefunction h in the presence 

of the Neumann condition at r = 1///. In the souce-probe calculation, the source graviton had 
Q oc rus and the probe graviton had mass m = nip. 

Moving to the massive ABJM field theory, fi is the same deformation parameter. The integer k 
can be used to define the 1/R scale with respect to some other scale (as, for instance, in the fuzzy 
sphere vacuum where R = Rgph/k, with Rsph the fuzzy sphere radius). As we noted before, there 
is no Mp inherent to ABJM since pure ABJM is of course conformal but Mp does appear as a 
unit of length in the same way as in AdS/CFT. For example, the physical radius of the spacetime 
sphere corresponding to the worldvolume fuzzy sphere is i?^^^ = ^^^Tr (^X^xj^ = 2Trkfj,Nlp. 
Then mg and rUp, or rather Ng and Np, should be obtained as sizes of the ABJM solutions 
corresponding to gravitons, i.e. sizes of vortices, with Ng + Np < N. Unfortunately, for 
technical reasons we had to restrict ourselves to the case of = 2 but we should, in principle, 
consider the case 1 <^ Ng, Np <^ N. It is not clear if it is sufficient to consider A^^ vortices on 
top of each other (i.e. an charge Ng vortex) , or if we need something more sophisticated in the 
case 1 <^ Ng <^ N. Finally, we should have a corresponding object for mi and the graviton 

wavefunction h should, in principle, match on to a vortex wavefunction. One possibility for 

the latter would be to define it as a(r) — n or g{r) — g{oo) = g{r) — 1, which decays to zero with 
the mass of the smallest excitation of the theory, i.e. as e~^^^ in the case of the A^^ = 1 vortex. 
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• Vortex calculations: The classical scattering of vortices for parallel separation gave 

1 

■j^Laass = ^{iJ- + /^y ) - XI nKQ{2nrij)vl^i ij , (8.7) 

i ij 

for N = 2, i.e. in the U{2) x U{2) ABJM theory. If we wanted to parallel the computation of 
BFSS, the next term to calculate would be at one-loop. Specifically, to obtain the interaction 
potential we must consider the classical solution of two vortices, one at 2; = and one at 
= 6* + v^t, and then calculate the one-loop determinants Z^„Ci;„/2 around it. To obtain 
matching with our tree-level graviton scattering result, we would need at least that the one- 
loop determinants give a result proportional to msiUpv'^, but even that seems very difficult to 
check directly, and we leave it for further work. 

In the case of transverse vortex separation, we anticipate the same situation as in BFSS in 
that, if we take two vortices in different SU(N) blocks then at the classical level they are non- 
interacting, even if they have non- vanishing velocities. Even though we will run into the problem 
of regulating their infinite energies, it is instructive to see how this works in the pure ABJM 
case. There, for a single vortex, = v^I for A = 2,3, 4. Therefore, by setting 

Y^{t) = {xf + vtt) liVixJVi + {xi + v^t) 1n,xN, , (8.8) 

and the vortices at the same worldvolume positions in the ABJM action, we obtain the nonrel- 
ativistic mass of the supergravitons exactly as in BFSS, 

,2 



Eo + {Ni + N2)^, (8.9) 



where Eq is the energy of the free vortices at rest and the interaction potential comes from the 
1-loop fluctuation around the vortex solution with X^{t) = U + v'^t. There are two immediate 
subtleties that arise: Firstly, it is not clear whether this will work for both X* being the four 
fuzzy sphere directions, as well as for the four remaining transverse directions and, secondly we 
have no explicit realization of these vortices in the massive ABJM model. Both of these are left 
for future work. 

Gravity calculations: Even if we could match our graviton scattering result (single graviton 
exchange) 

Knt = rUp^^h^^ OC msnipV^^i^^p ^^reZ.sp > (8-10) 

with a corresponding vortex computation like the one-loop one suggested above, it remains to 
find a gravity calculation matching the classical vortex scattering potential. For that, we need 
a calculation that will give a result proportional to vf^i -j . It does not necessarily need to be 
a classical calculation but it should certainly be something different from the leading single- 
graviton exchange calculated above a la 't Hooft, i.e. we want to find some term unaccounted 
for by that calculation. 

One possibility suggests itself immediately. Since the graviton itself should correspond to a 
vortex, it is not unreasonable to expect that we should also be able to find an effective 3- 
dimensional Lagrangian for interacting supergravitons. Its construction, however, is plagued 
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with subtleties. The supergraviton (a DO-brane) and the other £)-branes of the theory are 
also solutions of the 10-dimensional supergravity action, but only after introducing explicit 
delta function source terms in it. So we should, in principle, find a solution for interacting 
supergravitons and plug it back in the action, exactly as we did for the corresponding vortices. 
But that of course means that we now also need to do something about the extra dimensions. 
The most conservative route would be to explicitly perform the integration over the extra 
dimensions. Finally, one must still subtract the nonzero result for the background without 
supergravitons, corresponding to the ABJM field theory without vortices. To this end we start 
with 

^11 = ^EH + ^matter + ^source — {^EH + ^matter)bgr ■ (^-H) 

To arrive at the effective Lagrangian, we evaluate Cu on its supergraviton solution, and integrate 
out over all the directions (parallel and transverse), obtaining 

Leffmt)) = J (fx(fy^Cu{h^4y,r;yiit))), (8-12) 

where 

9iJ.u = g^^J + hfj,^{x, y; yi{t))6{xu - t) , (8.13) 

and wc drop the S{xii—t) factor from the metric together with the integration over xn, following 

[15]. Note that the 5{x~) factor in = h S{x~) can be eliminated from ^e// in this way 

only if it is the same for all the gravitons under consideration. But the different velocities of 
different gravitons mean that they do not move in exactly the same direction, and the delta 
function cannot be eliminated! So this approximation is exact only for the first order in the 
expansion in v, i.e. for the v'^ term we are interested in, but there will be extra terms starting 
at order v^. In order to calculate such effects, one must consider the full interacting problem, 
with differently oriented Shockwaves, 't Hooft showed that in flat space this is reproduced by 
the source-probe calculation used by [15], but now in our curved background the calculation 
might be incomplete. Nevertheless, this proposal passes several simple tests: 

— On the background (in the absence of a supergraviton) we have, correctly >Cii=0, since 
jC-source = by definition. 

— In the presence of a single (free) supergraviton, Ceh — {C'EH)bgr is also zero, since the 

only nonzero equation of motion is G — (G )i,gr oc A/i , and g = 0. The only 

contribution to L ejj comes from the source term, the integral of Csmircei 

which is exactly 

equal to the energy of the free supergraviton, as it should be. 

— In the case of interacting supergravitons however, we will have a nonzero interaction energy 
coming ixoTo. llEH—i^C.EH)hgr- Again, this energy will only be nonzero for a relative velocity. 
If the gravitons are parallel, it will be as if we have a single graviton of source equal to 
the sum of the all sources, for which we get no interaction. Further, Galilean invariance 
imples that the interaction has to enter as integer powers of (gj — qj)'^ = v^^i j^, so that we 
will only get v^, u^, f^, ... terms. 

^®Here x~ is DLCQ time, but we rewrite in terms of usual time t = xq and trivially integrating over xu at fixed 11th 
momentum of the fields. 
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— And finally the 1-loop ABJM interaction appears from classical corrections to the 6{x~)- 
factored Lagrangian as it should (although, of course, there must also be quantum gravity 
corrections as well, i.e. one must consider in the general the quantum effective Lagrangian 
instead of the classical Lagrangian). 



8.3 M5 brane backreaction 

As a final point in this section, we now sketch how to construct the solution for adding the backreacted 
M5-brane correction to the background corresponding to massive ABJM. As we saw in Section 5, in 
the type IIB theory the background corresponds to a pp-wave, and the M5-brane corresponds to a 
D3-brane wrapping the defined by ri = fixed, r2 = inside the maximally supersymmetric pp- 
wave. Such a solution has been implicitly written in Fig. If of LLM [42]. Moreover, if the M-theory 
dual giant graviton is to have size R = l/n, so too must the D3-giant in IIB. Our strategy will 
therefore be to (a) use the LLM prescription to find an explicit solution, (b) T-dualize this to type 
IIA and (c) lift to M-theory. Recall that the LLM solution for any i BPS type IIB supergravity 

ds'^ = -h-^{dt + Vidx'f + h'^{dy'^ +dx'dx') + ye^dnl + ye-^dnl 

h"^ = 2?/ cosh G 

ydyVi = eijdjZ 

y{diVj - djVi) = eijdyZ 

z = ^tanhG (8.14) 

F = dBtA{dt + V) + BtdV + +dB 

F = dBtA{dt + V) + BtdV + dB 

Bt = -^yV«; Bt = -\y'e-'^ 

dB = ~^y*sdC-^y, dB = -]y'ndi'-^) 
4 4 

is completely determined by a single function z which is itself the solution of the linear equation 

dAz + ydy(^^^ =0. (8.15) 

and absence of singularities means that z{xi,X2,y = 0) = ±1/2, and z ■H- —z is particle-hole duality. 
Note that the equation for $ = z/y'^ is the Laplace equation in 6d flat space, with 4d spherical 
symmetry, with y the radial coordinate for these 4 dimensions. The solution is 



zixi,X2,y) = -TT t^/ n- ' 1 2 + ^ 
27r Joj, [x - x'Y + y^ 

Vi{xi,X2,y) -■ 



(8.16) 



dx'^ 



27r/a^, (f-x')2 + : 



where a is a contribution from infinity if z is constant outside of a circle of very large radius (asymp- 
totically AdS^ X 5*5 geometries), with a = ±1/2 if z = ±1/2 asymptotically, P is a droplet, and rii is 
the unit normal vector on it, pointing towards the region of z = +1/2. 
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The solution we want is a superposition of two simple solutions: 

l)The pp-wave: This solution is obtained by having the boundary condition at y = as a 
half-filled plane, i.e. z{x'^,X2,Q) = ^sga{x2). Explicitly performing the integrals above gives 

1 X2 

(8.17) 



z{x2,y) = --===, 
2 ^/x^ + y2 



1 1 

V^ = ■ ^9 = 

To put the geometry into the standard pp-wave form 

ds"^ = -2dtdxi - [rl + rl)dt^ + dfl + drf , (8.18) 

we defined ^ 

y = rir2; X2 = -{rl-rl). (8.19) 

2) The AdSr, x solution: Now the boundary condition at y = is a large spherical droplet, i.e. 
for 5 = z — 1/2, the 6d Laplace equation for $ = z/y^ has source on a disk of radius tq. Integrating, 
one obtains 

~t . r'^-rl + y'^ 1 



r sin (j)Vi + r cos (j)V2 = V^{r, y; ro) 



9 I 9 I 9 
I- + .(/- + /•„ 



Defining new coordinates by 



2 \\/{r^ + + y'^Y - Ar'^rl 



y = vq sinh p sin 6 

r = ro cosh p cos ^ (8-21) 

d) = 6-t 



puts the metric into the more familiar form, 

ds'^ = ro - cosh^ pdt'^ + dp^ + sinh^ pdO| + dO'^ + cos^ 6*^^^ + sin^ OdCll , (8.22) 

with the area of the droplet in the {xi,X2) plane A = 47r^/p and where Ip = g^l'^^fd . This is clearly 
proportional to A^, the number of branes. Large N therefore corresponds to a large droplet. That 
the pp-wave is the Penrose limit of AdS^ x translates into the fact that near the surface of a large 
spherical droplet, locally the droplet looks like a half-filled plane. 

Since the equation for the only unknown function z{xi,y) is linear, we can trivially construct a 
solution that on the boundary give superposition of circles (and maybe half-filled planes), by just 
adding up the solutions for z{xi,y). In particular, for concentric circles we get 

z = J2i-^y^'~^{r,y;rS^) > 

i (8.23) 

^0 = E(-ir"''^^(-'y;^oO • 
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To obtain the solution we want - a D3-brane wrapping the with ri = ri^o, r2 = and {t,Xi) 
fixed inside the pp-wave - we translate into LLM language which, from (8.19), means that we need 
X2 = ^^^2^ ; y = and {t,xi) fixed. This is a small circle at a position X2 = away from the 

X2 = boundary of the half-filled plane. Specifically, we are interested in the circle of size vi^q = 
(or, in the notation used here, with n = 1, ri.o = 1). Thus the solution we seek is the (finite energy) 
solution of Fig.Sf in LLM, a distance X2 = 1/2 from the Fermi surface, and size given by how many 
branes we want wrapped, in this case one. 

In terms of the function z{xi,X2,y) this is just the sum of the pp-wave solution (8.18) and a 
rescaled version of (8.21), of area given by A?^ = 1 (one M5-brane), where we replace r of (8.21) with 
j,T — xq\, and xq = (0, ^). Unfortunately, this solution looks rather complicated, but to obtain the 
backreacted M5-brane we still need to T-dualize it to type IIA and lift up to M-theory. This, we 
leave as an interesting open problem. 

9 Conclusions 

In this paper we have looked for a Matrix theory-type construction from ABJM model. We have 
identified a BPS vortex solution of ABJM as corresponding to a spacetime supergraviton (a point- 
like object carrying DO-brane charge). In the case of the pure ABJM model, the identification is 
problematic since, in addition to the supergraviton, we also have also an infinite energy 2-brane to 
worry about. We find that this is resolved in the case of the massive maximally supersymmetric 
deformation of ABJM since, here, the solution has finite energy. We have also found a BPS kink 
solution of ABJM, using a certain limiting procedure. This solution can be identified with a Z)2-brane 
with one direction parallel to the ABJM worldvolume, and one direction transverse. The vacuum 
supersymmetric fuzzy sphere solution had already been identified with the D4-brane wrapped on S'^. 

In the latter half of the article, we identified the spacetime background corresponding to the 
massive ABJM deformation, and computed the leading graviton scattering interaction potential using 
a shockwave-probe method previously developed for the BFSS Matrix theory. We then computed the 
classical interaction potential of two vortices in the N = 2 ABJM model, and found that it differs from 
the supergraviton potential calculated by scattering two Shockwaves a la 't Hooft. We then speculated 
on how one could fix this mismatch, noting that there are many differences with respect to the BFSS 
computation. Chief among these is the fact that now, the two spatial worldvolume directions are 
special. We argued that perhaps one needs to find an effective Lagrangian for supergravitons that 
would match with the effective Lagrangian for vortex interaction, and that our gravity calculation 
could be matched by a one-loop calculation in the vortex background, similar to the BFSS case. 
Ultimately, though, it is clear that much more work is needed to find a concrete Matrix theory 
construction for ABJM. 
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A Graviton wavefunction equation 

In this Appendix we present some of the details of the solution of the graviton (shockwave) wave- 
function equation, 

A4>__ =Q<^^(x,y). (A.l) 

Wc begin by splitting the problem into two cases depending on whether the 11th dimension is par- 
allel to the worldvolume (transverse speparation) or perpendicular to the worldvolume directions 
(longitudinal separation). 

• Transverse separation: 

In this case, _L denotes directions transverse to the wave. Among the y = t,x,y directions, 
however, only y has this property, so eq.(A.l) becomes 

iJ-V3(f) + dl] {y, x) = QS{y)S\x) . (A.2) 

After a Fourier transform in the worldvolume direction y, 

h — {x,y) = dp —(pp{x), (A.3) 

we get the equation 

(1 - /x^f^)-^/' [dl - /(I - ^^(f) = QS^S) . (A.4) 

• Parallel separation: 

Now, there are two directions, y = (x, y) that are transverse to the wave and parallel to the 
3-dimensional worldvolume. Expanding the Laplacian and simplifying gives 

(1 - /X2f2)-V3 [^2 _ ^2(1 _ ^2^2)] ^^^^ ^ g^7(^) (^ 5) 

where now 



(x, y) = I d'p ^(t^p{S) . (A.6) 

We can actually treat both cases in a unified way since, in either case, we can write the equation that 
needs to be solved as 

" d 



, 2 2->2 2 

-p 



^p{x) = -Q5\x) . (A.7) 



Two points to note about this equation are (a) that it is manifestly spherically symmetric and (6) 
that if we put (5 = and p = ip we get the isotropic harmonic oscillator in d dimensions. To obtain 
a solution, we need to first impose boundary conditions, either at a finite distance or at infinity. The 
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natural one would be to impose normalizability of the wave function at infinity, but as we will see, 
this is not a good condition. 

To proceed, let's first solve the equation for Q = 0, using techniques for central interaction 
potential. Specifically, expand (pp in spherical harmonics as 



where the n-sphere spherical harmonics Y^""^ satisfy. 
Consequently 

f_ _ k{k + n-l)+n{n-2)/4 ,90 _ o 

TO 9 'T' H' P ' P 



yp,k{r) = . 



(A.8) 



(A.9) 



(A.10) 



Now note that when p = k = Owe obtain r"/^ and "/^ as solutions, since then <pp{x) has a constant 
and r^^" as solutions. Next, on setting p = ip, eq.(A.lO) takes the form 



y" nV- br^y + cy = Q. 



(A.11) 



This we recognize as a confluent hypergeometric equation whose solutions are, of course, (a linear 
combination of) confluent hypergeometric functions iFi and 2F1, 



y(r) = e 2 r"^ 2 



1 vm^ 



(A.12) 



On substituting a = ^1+"*^" ^) ^ 5 = ^2^2 j ^ _ ^2 ^j^^^ 



_ /c n + l, n + 1 ^2 



+C2 2F1 



p k n + l n+l 2 



(A.13) 



Since we have a spherically symmetric source and we are looking for a solution that goes to zero at 
infinity, we need to choose a spherically symmetric solution, i.e. k = 0. Also replacing p = —ip, we 
finally get 



(/>p(f) = e" 



4/i 4 



ip n + 1 n + 1 



; -ifipr^ ) +C2 2Fi{^ + 



ip n + 1 n + 1 



-; —ifipr 



Further, using the relation between 2-^1 and iFi 



2Fi{a,P;z) 



r(l-/3) 



T{l + a-P) 
we find that near r = 



iFi{a,P;z) + 



r(i-/3) u- 

r(a) 



iFi(l + a-/3,2-;3;z) 



(A.14) 



(A.15) 



3 — n3 — n 2 
47 + ^'^'"'^^' 



(A.16) 
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with the constant fixed as a function of the source Q (up to now neglected) as fohows. Near r = 0, 
the fuh equation, with nonzero source Q becomes 



which is solved by 



dlM^) = -Q6^{x) 

Q 



(j)p{x) = 



and so the dominating full solution near r = must be 



<t>p{x) = 



Q 



ip 3 — n 3 — n 



—ijipr 



(A.17) 
(A.18) 

(A.19) 



Notice that this solution behaves like l/r** ^. To this we can add the subleading correction (which 
behaves like a constant at near zero), 



^„.(f) = K,(-w)^eW'-V2 ( ^ + !i±l. 



(A.20) 



where Kp is an arbitrary p— dependent constant and we have also kept an explicit power of p above 
in order to emphasize the z = —ijipr^ dependence. 

Now to impose boundary conditions. Wc first try normalizability at infinity and find that both 
solutions behave in the same way as r — )■ oo, namely 



QT \ e-'^{-iij,pr'^)^'t^{-ifip)~^ 



(j)p{x) 



(n - l)^n-l 



{4(1^ A ) 



e 2 [inpr ) 4m {iixp) 



3 - II 
4 



ip I 3— n 



4m 



+ 



OC 



n + 1 n+\ 1 

T 2 r 2 



<t>p,2{x) 



e 2 



ip n + 1 n + 1 



4/i ^ 4 



2 ' 



—ijipr 



(...)r2M 



n+l 
2 



+(...)e 2 r 2^ 2 OC 



n+l 
f 2 



and correctly going to zero at infinity. But this behaviour is non-normalizable, since 



2 



dr 



Ini?. 



(A.21) 



(A.22) 



We can extract more information from the solution by taking the p — )■ oo limit of the result and using 
the fact that the confluent hypergeometric function 



COS 



2\/^-7r(^ - ^) 



as |a| oo. This, in turn means that as p — )■ oo 

4>p{x) C:^ 



■—^{±ipr) 2 COS 



^/^^{n - l)^n-i r'^' 



lb ipr + TT- 



n - 2 



(A.23) 



(A.24) 
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To continue, we need to fix Kp by imposing appropriate boundary conditions. Again, we distinguish 
two cases. 

Case 1. We first assume that we can put Kp = 0, thus ignoring <f)p^2{x)- Then, in the transverse 
separation case where n = 7, we obtain the graviton wavefunction at y = by a simple integration 
of ^p(x) over p. Assuming also that the large p region dominates the integral, we obtain 

^ r^ (A.25) 



C = lim ^, Z/^ (±i) 2 / dzz 2 cos 

-1 7-( 



n^7 ■\/7r(n — l)^r, 



n-2 
±iz + TT — : — 



which matches the fiat space case (corresponding to the BFSS analysis). Also note that there are no 
other dimensional parameters in the Q/r'' behaviour, it is just multiplied by a number (part of the 
number is an integral over the variable z = pr), for which we must take a limit since r(— 2) = 00. 
This independence of the result from /x is most likely the result of an incorrect initial assumption 
about the boundary conditions. 

In the case of parallel separation, n = 6 and the Fourier transform is 

/j2 poo r27r jg noo 

^e'^ct>p{x) = pdpj^ —e'vy---'ct>p{x). = j^ p dp Mpy)M^) (A.26) 

Here, by contrast, we are interested in the r — t- limit. If however, p stays finite, we find that the 
solution (A. 19) becomes p— independent, and then the above Fourier transform gives zero (or rather, 
S{y)). So we must again consider the large p limit in (A. 24). This region of large p can still give a 
nontrivial contribution to the p integral, since pr is not necessarily small. Then, also substituting 
n = 6, gives 

(j){y, r ^ 0) ~ -TTt) — 3 / '^P MPy) cos{ipr + vr) = -zr^ — 3 / dp p^ Jo{py) cosh(pr) 
LMiQr Jq Loiifjr Jq 

- -^-A^F^ (2,5/2; . (A.27) 



To understand what the behaviour of the solution near r = is, the identities 

2F,{a,f5;r,z) = w^lw^ ~ (-^)""2Fi(«, a + 1 - r,a + 1 - P;-) (A.28) 

r(p)r(7-a) z 

can be used to show that as y/r — > 00, the hypergeometric function 2-F1 ^2, 5; 1; — ~ — (^^nr Si-nd, 

subsequently ci>(y, r — >■ 0) ~ -^-^1 — r~rr7T iiear r = 0. Note then that we can actually set r = 

bile r' {y/r)^^ 

directly and obtain ^(r = 0) = 0, which is clearly not what we wanted, so again the initial boundary 
condition was incorrect. 
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Case 2. We now look for a more physical boundary condition, and we concentrate on the parallel 
separation case. In reality, the r space terminates at r = 1/fJ,, so we must impose a boundary 
condition there. This is an apparent singularity that one should be able to continue through, thus 
the appropriate boundary condition at r = is Neumann, 4''p{r = 1/ fi) = 0. The reason is that 
then the point r = acts as the origin in angular coordinates, for which the above condition is the 
only one that makes sense (see for instance [43]). Thus imposing 



-^Mr)\r=l/^, = 0, (A.29) 



on the function 



(A.30) 



and using the identity 



we obtain the condition 



-^iFi{a,b;z) = ^iFi(a + 1,6 + 1; (A.31) 



1-^1 \4:tx ^ 4' 2' /i J ^ 7 \4tj. ^ Aj ^ 4 ' 2' /i J 

We already saw that at r = the Fourier transform of the first term in (A.30) vanishes, so we need 
to analyze only the second term. Then 

cj>{y, r = 0) = J ^e'P-yKp{-i,,pf/^ (A.33) 

and the right hand side of (A. 32) gives in the large p limit 

^^^ + 'f^PV' (A.34) 



16 1 - i V2/i, 

which means that we could close the contour of integration dp with a semicircle in the upper 
half plane (since Kp{—ifj,p)^/'^ goes to zero as \p\ — t- oo), and thus if the integral above would be 
one-dimensional instead of two, it would be given by the residues at the poles in the upper half p 
plane. Kp has poles in the upper half plane of the complex p, given by the solutions to the equation 

„ f ip ,7 7 -ip\ _ 4 Zip 7\ p «P , 11 9 -ip\ 

We can easily see, however, that p = is a solution to this equation, so half the residue at p = will 
contribute to the integral, giving a constant contribution (independent of y) instead of an exponential. 

The next solutions that we obtain for the poles (with Mathematica) are p = 0. i: 9.14066z/x, for 
which the residue does indeed give an exponential. So the main contribution to the integral gives a 
constant, plus an exponential ~ e~"^i^. 

Here mi the lowest imaginary part in the upper half plane among these solutions. Note that 
we cannot say for certain that mi = 9.14066/i, since Mathematica searches are always in a neigh- 
bourhood; if there is a solution of small imaginary part at large real part, we cannot say from the 
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numerical result. Continuing the numerical search for solutions, p, of the root equation, we find 
0. ± 13.4767z/x; 0. ± 17.6461z/x; 0. it 21.7516z/x, with a rather general starting point (even a real one). 
This seems to hint that all solutions are pure imaginary? If so, we would have a chance to prove that 
mi = 9.14066;u. 

Thus if the p integral were be one-dimensional, we would get (f){y — t- oo,r = 0) ~ e"™'^^. We can 
then at least conclude that if we restrict the dependence to only one of the y coordinates, we indeed 
get 4>(y — )■ oo,r = 0) ~ e^™^^. A similar result is expected for a two-dimensional y, since then we 
would need to do the integral 

/•OD 

Hlvlr = 0) = / pdp Jo{py)kj,{-ipipf^ , (A.36) 
which we cannot perform but the small y behaviour would likely be given by expanding Kp{—ifj,p)^/'^ 

POD 

at large p, which gives oc / dp p~^Jo{py) oc y'^. It is not very clear how to obtain the large y 

^0 

behaviour, especially since as we saw above, we are interested in a subleading term. The leading term 
is indeed (rather easily) obtained as a constant, since \p\ — > 0, 

and 

dv 

<^(|y|,r = 0)oc / ^Mpy) = const. (A.38) 
Jo P 

(the exact form of the constant we cannot be sure of, since the integral above is outside the range 
of validity for formulas we could find). It is very likely then that the same massive behaviour e~^^y 
persists for the first subleading term but this needs to be checked. 



References 



[1] J. Bagger and N. Lambert, "Modeling multiple M2's," Phys. Rev. D75 (2007) 045020, 
arXiv : hep-th/061 1 108. 

[2] J. Bagger and N. Lambert, "Comments On Multiple M2-branes," JHEP 02 (2008) 105, 
arXiv : 07 1 2 . 3738 [hep- th] . 

[3] J. Bagger and N. Lambert, "Gauge Symmetry and Supersymmetry of Multiple M2-Branes," 
Phys. Rev. D77 (2008) 065008, arXiv : 0711 . 0955 [hep-th] . 

[4] A. Gustavsson, "Algebraic structures on parallel M2-branes," Nucl. Phys. B811 (2009) 66-76, 
arXiv: 0709. 1260 [hep-th]. 

[5] O. Aharony, O. Bergman, D. L. Jafferis, and J. Maldacena, "N=6 superconformal 

Chern-Simons-mattcr theories, M2-branes and their gravity duals," JHEP 10 (2008) 091, 
arXiv: 0806. 1218 [hep-th]. 

[6] H. Nastase, C. Papageorgakis, and S. Ramgoolam, "The fuzzy S'^ structure of M2-M5 systems 
in ABJM membrane theories," JHEP 05 (2009) 123, arXiv : 0903 . 3966 [hep-th] . 



38 



[7] H. Nastase and C. Papageorgakis, "Fuzzy Killing Spinors and Supersymmetric D4 action on the 
Fuzzy 2-sphere from the ABJM Model," JHEP 12 (2009) 049, arXiv: 0908. 3263 [hep-th] . 

[8] T. Banks, W. Fischler, S. H. Shenker, and L. Susskind, "M theory as a matrix model: A 
conjecture," Phys. Rev. D55 (1997) 5112-5128, arXiv:hep-th/9610043. 

[9] J. Castelino, S. Lee, and W. Taylor, "Longitudinal 5-branes as 4-spheres in matrix theory," 
Nucl. Phys. B526 (1998) 334-350, arXiv:hep-th/9712105. 

[10] H. Grosse, C. Klimcik, and P. Presnajder, "Finite quantum field theory in noncommutative 
geometry," Commun. Math. Phys. 180 (1996) 429-438, arXiv:hep-th/9602115. 

[11] D. N. Kabat and W. Taylor, "Spherical membranes in matrix theory," Adv. Theor. Math. 
Phys. 2 (1998) 181-206, arXiv:hep-th/9711078. 

[12] J. Gomis, D. Rodriguez-Gomez, M. Van Raamsdonk, and H. Verlinde, "A Massive Study of 
M2-brane Proposals," JHEP 09 (2008) 113, arXiv: 0807 . 1074 [hep-th]. 

[13] A. Sen, "DO branes on T(n) and matrix theory," Adv. Theor. Math. Phys. 2 (1998) 51-59, 
arXiv : hep-th/9709220. 

[14] N. Seiberg, "Why is the matrix model correct?," Phys. Rev. Lett. 79 (1997) 3577-3580, 
arXiv : hep-th/9710009. 

[15] K. Becker and M. Becker, "On graviton scattering amplitudes in M-theory," Phys. Rev. D57 
(1998) 6464-6470, arXiv:hep-th/9712238. 

[16] G. 't Hooft, "Graviton Dominance in Ultrahigh-Energy Scattering," Phys. Lett. B198 (1987) 
61-63. 

[17] K. Kang and H. Nastase, "Planckian scattering effects and black hole production in low M(P1) 
scenarios," Phys. Rev. D71 (2005) 124035, arXiv:hep-th/0409099. 

[18] K. Kang and H. Nastase, "High energy QCD from Planckian scattering in AdS and the 
Froissart bound," Phys. Rev. D72 (2005) 106003, arXiv:hep-th/0410173. 

[19] H. Nastase, M. A. Stephanov, P. van Nieuwenhuizen, and A. Rebhan, "Topological boundary 
conditions, the BPS bound, and elimination of ambiguities in the quantum mass of solitons," 
Nucl. Phys. B542 (1999) 471-514, arXiv:hep-th/9802074. 

[20] D. E. Berenstein, J. M. Maldacena, and H. S. Nastase, "Strings in flat space and pp waves from 
N = 4 super Yang Mills," JHEP QA (2002) 013, arXiv:hep-th/0202021. 

[21] D. A. Lowe, H. Nastase, and S. Ramgoolam, "Massive IIA string theory and matrix theory 
compactification," Nucl. Phys. B667 (2003) 55-89, arXiv:hep-th/0303173. 

[22] J. Gomis, D. Rodriguez-Gomez, M. Van Raamsdonk, and H. Verlinde, "Supersymmetric 

Yang-Mills Theory From Lorentzian Three- Algebras," JHEP 08 (2008) 094, arXiv: 0806. 0738 
[hep-th] . 



39 



[23] S. Mukhi and C. Papageorgakis, "M2 to D2," J HEP Oh (2008) 085, arXiv: 0803. 3218 
[hep-th] . 

[24] M. Arai, C. Montonen, and S. Sasaki, "Vortices, Q-balls and Domain Walls on Dielectric M2- 
branes," arXiv: 0812. 4437 [hep-th]. 

[25] C. Kim, Y. Kim, O.-K. Kwon, and H. Nakajima, "Vortex-type Half-BPS Solitons in ABJM 
Theory," Phys. Rev. D80 (2009) 045013, arXiv: 0905. 1759 [hep-th]. 

[26] R. Auzzi and S. Prem Kumar, "Non-Abelian Vortices at Weak and Strong Coupling in Mass 
Deformed ABJM Theory," arXiv : 0906 . 2366 [hep-th] . 

[27] C. G. Callan and J. M. Maldacena, "Branc dynamics from the Born-Infeld action," Nucl. Phys. 
B513 (1998) 198-212, arXiv:hep-th/9708147. 

[28] P. L. H. Cook, R. de Mello Koch, and J. Murugan, "Non-Abelian Blonic brane intersections," 
Phys. Rev. D68 (2003) 126007, arXiv:hep-th/0306250. 

[29] G. W. Gibbons, "Born-Infeld particles and Dirichlet p-branes," Nucl. Phys. B514 (1998) 
603-639, arXiv:hep-th/9709027. 

[30] J. P. Gauntlett, J. Gomis, and P. K. Townsend, "BPS bounds for worldvolume branes," JHEP 
01 (1998) 003, arXiv:hep-th/9711205. 

[31] T. Fujimori, K. Iwasaki, Y. Kobayashi, and S. Sasaki, "in preparation,". 

[32] J. Gomis, A. .J. Salim, and F. Passcrini, "Matrix Theory of Type IIB Plane Wave from 
Membranes," JHEP 08 (2008) 002, arXiv: 0804. 2186 [hep-th]. 

[33] K. Hosomichi, K.-M. Lcc, S. Lee, S. Lee, and J. Park, "N=5,6 Superconformal Chern-Simons 
Theories and M2-branes on Orbifolds," JHEP 09 (2008) 002, arXiv: 0806. 4977 [hep-th]. 

[34] T. H. Buscher, "Quantum corrections and extended supersymmetry in new sigma models," 
Phys. Lett. B159 (1985) 127. 

[35] T. H. Buscher, "A Symmetry of the String Background Field Equations," Phys. Lett. B194 
(1987) 59. 

[36] T. H. Buscher, "Path Integral Derivation of Quantum Duality in Nonlinear Sigma Models," 
Phys. Lett. B201 (1988) 466. 

[37] S. Ramgoolam, "On spherical harmonics for fuzzy spheres in diverse dimensions," Nucl. Phys. 
B610 (2001) 461-488, arXiv:hep-th/0105006. 

[38] B. Janssen, Y. Lozano, and D. Rodriguez-Gomez, "Giant gravitons in AdS(3) x S**3 x T**4 as 
fuzzy cylinders," Nucl. Phys. B711 (2005) 392-406, arXiv:hep-th/0406148. 

[39] K. Hanaki and H. Lin, "M2-M5 Systems in N=6 Chern-Simons Theory," JHEP 09 (2008) 067, 
arXiv : 0807 . 2074 [hep-th] . 



40 



[40] A. Jevicki and H. Nastase, "Towards S matrices on flat space and pp waves from SYM," 
arXiv : hep-th/0501013. 

[41] K. Becker, M. Becker, J. Polchinski, and A. A. Tseytlin, "Higher order graviton scattering in 
M(atrix) theory," Phys. Rev. D56 (1997) 3174-3178, arXiv:hep-th/9706072. 

[42] H. Lin, O. Lunin, and J. M. Maldacena, "Bubbhng AdS space and 1/2 BPS geometries," JHEP 
10 (2004) 025, arXiv:hep-th/0409174. 

[43] E. Witten, "Anti-de Sitter space, thermal phase transition, and confinement in gauge theories," 
Adv. Theor. Math. Phys. 2 (1998) 505-532, arXiv:hep-th/9803131. 



41 



